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We demonstrate multiferroic behavior in trimerized Mott insulators through the interplay between spins and 
electric dipole moments resulting from electronic charge fluctuations in frustrated units. The model consists of 
stacked triangular layers of trimers with small intertrimer exchange interactions J' and J'\ Ferroelectric states 
coexist with ferro- or antiferromagnetic orderings depending on the value of the magnetic field H and the sign 
of the interlayer exchange J'\ The electric polarization undergoes abrupt changes as a function of H. 
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Introduction. — Frustrated Mott insulators have been the fo- 
cus of research during the last decades 1 1 1. The combination 
of geometric frustration with strong electron-electron interac- 
tions often leads to unusual collective behaviors. For instance, 
geometric frustration is a precondition for having nonuniform 
electronic charge distributions in the Mott phase of half-filled 
Hubbard models |2|. This requirement implies that the lattice 
must contain triangular units of magnetic ions. The simplest 
unit is an equilateral trimer of S = 1/2 ions, like Cu^^, whose 
ground states are two S = 1/2 doublets. This is true as long 
as the spin is a half-integer. Since the spin quantum number 
is not enough to characterize each ground state, it is necessary 
to introduce an efiTective orbital degree of freedom (DOF) de- 
scribed by a T = 1/2 pseudospin variable |2|. While (Tx, Ty) is 
proportional to the efiTective electric polarization operator, 
is proportional to the orbital magnetic moment produced by a 
current density that circulates around the trimer |2| . 

Materials that naturally include spin and orbital DOFs ex- 
hibit a variety of complex behaviors f3\. In particular, or- 
bital ordering often reduces magnetic frustration by creating 
disparities between efifective exchange constants (H. While 
the J-orbital DOF of transition metals is related to electronic 
charge distributions with diflTerent quadrupolar moments, the 
orbital DOF of triangular molecules carries a net electric 
dipole moment. Thus, collective behaviors of coupled trimers 
can lead to multiferroic phenomena or magneto-electric ef- 
fects arising from the interplay between the spin and orbital 
DOFs. 

Trimers of magnetic ions are rather common in organic and 
inorganic compounds II5HT31. They also exist in crystalline 
systems such as spin tubes, but the inter- and intratrimer ex- 
change interactions are of comparable magnitude 1 14-161. 
Although La4Cu3MoOi2|[I21 [TSl is an ideal realization of 
weakly coupled trimers, the trimer superlattice does not fa- 
vor a ferroelectric ordering |[T9ll. The advantage of organic 
environments is their flexibility for designing specific trimer 
lattices by choosing adequate ligand fields. While intramolec- 
ular exchange in frustrated molecules is a current topic of fo- 
cus in magnetochemistry, little is known about the collective 
behaviors induced by intermolecular couplings. 

The purpose of this Letter is to bridge the gap between 
molecular and crystal magnets by demonstrating that multi- 
ferroic collective phenomena can arise from intertrimer ex- 



change. After noticing that each trimer carries an internal 
electric dipole moment, it is natural to ask what is the effective 
coupling between these moments. The answer depends on the 
nature of the trimer lattice and the sign of the exchange inter- 
action |19, 20|. By demonstrating that a trimerized triangu- 
lar lattice leads to ferro- or ferrielectric ordering, we provide 
guiding principles for designing new multiferroic materials. 

We start by considering a Hubbard lattice of stacked tri- 
angular layers of trimers with small intertrimer hopping and 
large on-site Coulomb repulsion U. We find multiferroic 
ground states that remain stable up to magnetic fields above 
which the magnetization (M) and the electric polarization (P) 
change discontinuously. These multiferroic states and strong 
magneto-electric effects are direct consequences of the effec- 
tive interaction between spin and orbital DOFs. 

Model — The half-filled Hubbard model on the trimerized 
stacked triangular lattice of Fig.[T]is 

= - Z ^'V {^Uj. + H.C.) + I J] („, - 1)2 , (1) 

ijo- i 

where cj^ (c.^) is the creation (annihilation) operator of an 
electron with spin cr at a site /, ni = 2o- ^^L-^/o- number 
operator, and the hopping amplitudes are tij = t for / and j 
in the same trimer, and tij = f(f') when / and j are nearest- 
neighbor sites belonging to different trimers of the same layer 
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FIG. 1. (color online). Trimerized triangular lattice. Shaded triangles 
represent trimers with dominant hopping amplitude. The layers are 
stacked along the c direction. A,B,C are the indices for the three 
trimer sublattices. 



(adjacent layers). In what follows we consider the strong 
coupling limit U :$> t, though we will also comment on the 
intermediate-coupling regime U > t. 

When U ^ the half-filled Hubbard model can be reduced 
to a Heisenberg Hamiltonian *Hspin by applying degenerate 
perturbation theory to the second order in tif. 



*^spin — J ^ ^ ^ ^fi,r ' ^T],r ~^ J ^ ^ ^ ^ri,Y ' S^,r+e;, 

-h r ^ S;^,r • S;7,r+U3 - g/^BH ^ s'^^^. (2) 



r,fi>r] 



r,T] 



Here we have added a Zeeman term to include the eff'ect of 
an external field H (g is the gyromagnetic factor and yUg is 
the Bohr magneton). We have also refined our notation by 
introducing the trimer coordinate, r = Zi= 1,2,3 ^/U/, with U/ 
being primitive vectors for the trimer lattice and /d,r] = 1, 2, 3 
denoting the three ions of each trimer. are relative vectors 
between intralayer nearest-neighbor trimers (see Fig.[T] ). *Hspin 
also describes systems of coupled spin tubes [ 2TM26I . 

The reduction of to a spin Hamiltonian for » 1 sug- 
gests that only magnetic DOFs remain active at low energies. 
However, virtual charge fluctuations can produce electric cur- 
rents in loops or electric dipoles |2|. We will see below that 
this is indeed the case for the ground states of as long as 
t » The exchange constants /, and f are propor- 

tional to t^/U, f^/U, and f'^/U, respectively. Therefore, spin 
trimers are weakly coupled (/ » J\J'') for t » t\f\ and 
that will be the regime of interest from now on. Although *K 
leads to an antiferromagnetic (AFM) interlay er exchange 
we will also consider the FM case, which can be realized for 
superexchange paths through intermediate ions | 27 -29|. 

A single trimer has four degenerate ground states, namely, 
two 5 = 1/2 doublets, that can be labeled as \S^,t^)j., where 

and are the (±1/2) eigenvalues of 

2 

= + 4,r + 4,r' = y^^l'*- ^ ^2,r ' S3,r. (3) 

is the scalar spin chirality that is proportional to the efl'ec- 
tive current density operator in the trimer r. It closes an SU(2) 
algebra with the operators, 

= i [2S2,r • S3,r " Si,r • (S2,r + S3,r)] OC P^, (4) 



= — pSi,r • (S2,r " S3,r) OC P^, 

V3 



(5) 



which are proportional to the x and y components of the trimer 
electric dipole moment [2J . These spin and orbital trimer op- 
erators commute with each other: [r^, S^,] = 0. 

The S = 3/2 excited states of the single trimer can be pro- 
jected out as long as the single-trimer gap of 3// 2 is much 
larger than J' and \J''\. The low-energy efl'ective Hamilto- 
nian is simply *Keff = ^^^^inP, where P is the projector onto 
the subspace spanned by the direct product of single-trimer 
ground states |5^T^)r = 1+1/2, ±1/2)^ EDHH. Hff only in- 
cludes contributions from interdimer interactions, which leads 



to P^ii,Y ' ^r],r 'P = ^s^,r^ ' P^r],Y'P • H is then Convenient to in- 
troduce the projected spin operators: 



(6) 



where = (cos^^, sin^^,o) with = 4(jd - l)7r/3 (// = 
1, 2, 3) are parallel to the displacement vectors from the cen- 
ter of the trimer to each spin site, s^r must be proportional 
to Sr because this is the only vector under spin rotations in 
the single-trimer ground state subspace. The proportionality 
factor cannot include because it is odd under time reversal. 
This implies that *Heff will only include the and orbital 
variables associated with the local electric polarization of each 
trimer. When the spin operators of Eq. ([2]) are replaced by the 
projected operators of Eq. ([6]), we obtain 

r,T] 

X [l - 4Tr • 11;^ + 2Tr+e, * " 8Tr+e, ' nr^Tr ' Ur^] 
r 

-gfiBHYsi. (7) 



The continuous symmetry of the interlayer orbital coupling is 
a consequence of the C3 lattice symmetry. 

Above the saturation field. — We will first consider the case 
of magnetic fields that are large enough to polarize each trimer 
into the 5^ = 1/2 state (34]: H > H,^^. Since Sr • Sr+U3 = 
Sr • Sr+e^ = 1/4 holds in this case for the low-energy sector 
<Sfp of "Hqs, the projection Psfp of "Heff into <Sfp leads to 



r,r] 



const. 



(8) 



Here, the linear terms in r vanish because ^/u = 0- ^^eff 
is an eff'ective model for the orbital DOE which clearly shows 
that magnetic exchange interactions induce eff'ective exchange 
couplings between the electric dipole moments. According to 
Eq. ([S]), an AFM intralayer exchange, J' > 0, induces a fer- 
roelectric (EE) exchange between electric dipoles in the same 
layer. In contrast, a FM (AFM) interlayer coupling results in 
a EE [antiferroelectric (AFE)] coupling along the c direction. 
We will now consider both possibilities by using a semiclas- 
sical approach and assuming that magnetic and orbital order- 
ings are three- sublattice structures. Since the eff'ective inter- 
action is XF-like, we will assume that (rj) = 0, Vr. The XY 
components are determined by three variational parameters O/ 
with / = A,B,C being the sublattice index of the trimers (see 
Fig.[l]): (rj ) = (o-r/2) cos O/ and (r^) = (o-r/2) sin O/ for re/. 
Here, CTr = 1 [(Xr = (-1)"^ with being the layer index] for 
< > 0). The mean-field (ME) energy per site that 
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FIG. 2. (color online), (a) The MF phase diagram. The insets show 
M{H) and P{H) for J" < (left) and J" > (right). The dashed line 
shows an orbital transition due to the zero-point fluctuation. Also 
shown are schematic pictures of the MF solution [the arrows above 
(embedded on) trimers represent spins (orbitals)] : (b) the FM-FE, (c) 
FIM-FIE, (d) FM-AFEl, (e) FIM-AFE2, and (f) CAFM-FIE states. 
The solid (dashed) lines indicate in-layer antiparallel (parallel) or- 
bitals. Zero-field states are given up to an SU(2) global spin rotation. 



results from Eq. ([5]) is 

€h>h,,, = -t^ Yj cos (O/ - (D/+i) -h const. (9) 

l=A,B,C 

It is clear from this expression that the global minimum of 
^H>H,^, is obtained for O/ = O. This solution corresponds to 
a fully polarized FE (AFE) state for r < (r > 0); see 
Figs. [2jb) and (d). The arbitrary value of O implies that the 
MF energy is invariant under global orbital rotations along the 
z axis; i.e., there is U(l) orbital symmetry at the MF level. 

The next step is to introduce Holstein-Primakov bosons 
to describe the orbital fluctuations around the MF solution: 

+ ITr ' Cl' = V2t - Har Clr, " Cl' = al ^J2T - Har, and 

case for f > leads to the same "orbital-wave" Hamilto- 
nian given below by a trivial gauge transformation. In terms 
of the Fourier modes, a^^ and a]^, the Hamiltonian to order 1 /r 
is 



for f < 0, where Q = (cos O, sin 0, 0), 



(sin O, - cos 0, 0), Uar = a^a^, and r = 1/2. The 



^T,H>H,^ 



2 



k + 



6k = -(/-n-rk, 

Tk = ^ sin^ (O - (pr^) cos k • C;^ - y cos k^. (10) 



By performing a Bogoliubov transformation, we obtain 

'^rli>H,,. = 2:k[^k(4^k + 1/2) - ^k/2] with oj^ = ^el-yl. 
This dispersion relation has a zero mode at k = due to the 
U(l) invariance of the MF solution. However, after including 
corrections due to zero-point fluctuations, the energy density, 

= eH>H,,, + A6 with ^e = -(f - r')/3 + 2k ^k/2, is min- 
imized for O = (2n + l)7r/6, with n being an integer number. 
Quantum fluctuations reduce the number of ground states ac- 
cording to the sixfold symmetry of "H^^ and the maximum 
splitting between states with difl'erent values of O is of or- 
der 10"^ /^ Consequently, higher order corrections in our l/r 
expansion induce a very small gap in comparison to the ex- 
change constants. For this reason, we will say that the orbital 
k = mode is a quasi-Goldstone mode. 

States for H < Hsat- — Now we consider the general case 
of an arbitrary magnetic field. We propose a three- sublattice 
ordering in each layer and distinct even and odd layers for 
J'' > 0. We use (^p 0^) to represent (Sr) and assume (r^) = 
(1/2) cos Of, <T^) = (1/2) sin Of, and (rj) = for r g (/?,/) 
with b = e,odindl = A,B,C. 6'^ = 6^, = 0;^, and O^ = O;^ 
are assumed for J'' < 0. The MF energy density resulting 
from Eq. ^ is 

6mf = — ^ [sin ^f sin ^f cos (^f - 0f ) + cos ^f cos ^f ] 



b,l>l' 



[l-cos(<Df-<Df,)] 



+ — ^ [sin sin cos (^^ - 0^) + cos cos O'l^ 



[l.2cos(o^-0?)]-^X-^^/- 



(11) 



b,l 



The phase diagram obtained by minimizing ( pTj ) is shown in 
Fig. [2ja). As we already discussed, a FM-FE or FM-AFEl 
state is stabilized for H > ifsat depending on the sign of 
r'. When H < if sat, the MF ground state for f < is 
the ferrimagnetic (FIM)-ferrielectric (FIE) state of Fig. [2jc) 
with an in-plane "up-up-down" orbital ordering that accom- 
panies a corresponding up-up-down spin configuration, e.g., 
(OA,Ofi,Oc) = (0,0,0 + 7r)and (SA,Sfi,Sc) = (tti) (a 
permutation of sublattice indices does not change the energy). 
This state leads to 1/3 plateaus both in M{H) and P{H) that 
change discontinuously atH = //sat [see the inset of Fig.[2ja)]. 
The energy is minimized because the orbital ordering reduces 
the eflTective magnetic frustration. This is the reason why the 
ground state diflTers from the 120° structure for the spin-only 
AFM Heisenberg model on the triangular lattice. 

On the other hand, the zero-field solution for J'' > is the 
AFM-FIE state shown in Fig. [2jf) [the solution is given up 
to a global SU(2) spin rotation and 6i = 62 = 7i/2], which 
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is also a collinear spin state. In contrast to the f < case, 
a finite H induces the coplanar canted AFM (CAFM) state 
shown in Fig.[2jf), where M{H) increases Hnearly while P{H) 
remains constant at a 1/3 polarization plateau. Interestingly, 
the FIM-antiferroelectric (AFE2) state of Fig. |2je) is stabi- 
lized for f If < 4.635 in the intermediate field, leading to 
a 1/3 magnetization plateau as in the J'' < case. M{H) in- 
creases discontinuously while P{H) vanishes abruptly at the 
first-order transition to this state. The FIM-AFE2 state [see 
Fig. [2je)] is such that in-layer FE bonds are oriented along 
diff'erent directions on adjacent layers. This implies that there 
is an accidental extensive degeneracy at the MF level that is 
removed by fluctuations. However, the behavior of M(H) and 
P(H), namely, the large magneto-electric efl'ect, is the same 
for all the degenerated states. As for f' < 0, the mechanism 
for stabilizing these states is a suppression of magnetic frus- 
tration by orbital ordering. For any f > the transition to 
the FM-AFEl state [Fig.|2];d)] is of first order. 

Finally, we include the spin- wave analysis for H < if sat • 
We will mainly discuss the f < case (i.e., the FIM-FIE 
state). A more comprehensive analysis will be given else- 
where |35|. By choosing one of the FIM-FIE configurations, 
(Oa, O5, Oc) = (0, 0, 71) and (Sa, S^, Sc) = (T, T, i), we intro- 
duce the Holstein-Primakov bosons, Z?/r (air), for represent- 
ing the spin (orbital) fluctuations on the sublattice /. Spin 
and orbital fluctuations are decoupled in the linear spin-wave 
Hamiltonian. Moreover, the orbital Hamiltonian can be re- 
duced to Eq. ([To]) simply by redefining = air for each 
r G I = A,B,C. This observation implies the stability of the 
FIM-FIE state against orbital fluctuations, and the presence of 
the quasi-Goldstone orbital mode. The spin contribution of 
the spin- wave Hamiltonian is 

^5,FIM-FIE - 2 ZjV k'^-kMg* I \ # ) 

where = (^Ak, ^fik, ^ck) with /?/k being the Fourier trans- 
form of bir. P\i and 2k are given by 
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where = (f /9)Y.^[l ±2cos^ - (fir^)]e^^-^^ , 4 = 
(2/73) - r (1 - cos^,) s/ibH, and 4 = 4 + (2/73) - 
2giiBH. This Hamiltonian can be diagonalized by solving 
|(^k + 2k) (^k - 2k) -col^ = |36| and a positive-defined 
spectrum is obtained for any H < //sat, indicating that the 
FIM-FIE state is locally stable. The spectrum is gapped for 
< H < //sat. This also implies stability of the 1/3 mag- 
netization plateau for f > (i.e., the FIM-AFE2 state) at 
least for sufficiently large H and small f because both states 
are connected at f' = 0. Both spin and orbital fluctuations 
contribute to the zero-point energy that determines the precise 



structure of the orbital ordering. We find that the effect of spin 
ffuctuations competes against the effect of orbital fluctuations: 
the energy minima of the spin contribution are Omm = nn/S 
in contrast to Omin = (2n -h l)7r/6 for the orbital contribution. 
The ground state energy turns out to be almost insensitive to 
H/f but it is significantly affected by f /J\ which causes a 
phase transition from O^in = nn/S for the quasi-2D region 
-0.8 < r/J\< 0) to Ornin = (2n + 1)7t/6 for r/f < -0.8 
[see Fig.iga)]. 

Conclusions. — In summary, the interplay between spin and 
orbital DOFs of weakly coupled trimers leads to multifer- 
roic behavior and strong magneto-electric effects. Frustra- 
tion plays a fundamental role in different stages of this prob- 
lem. While it is crucial for the emergence of orbital DOFs 
that carry electric dipole moments |2|, it is also essential for 
stabilizing the multiferroic orderings depicted in Fig. [2] The 
quasi-Goldstone orbital mode of these ordered states can be 
indirectly observed in the low temperature (T) regime by mea- 
suring its (d is the spatial dimensionality) contribution to 
the specific heat for T higher than the tiny orbital gap. Spin 
excitations give a negligible contribution to the specific heat 
in the plateau phases because of the much larger spin gap. 

Our perturbative approach for weakly coupled trimers can 
be extended to the intermediate-coupling regime U > tSiS long 
as 1^1 » t\ The effective spin-orbital Hamiltonian forU>t 
must be derived directly from *H because "Hspin no longer re- 
produces the low-energy spectrum of 'H. The electric dipole 
of each trimer is much stronger in this regime (comparable 
to ea, where e is the electronic charge and a is the trimer lat- 
tice parameter) and can lead to polarizations of order lyuC/cm^ 
if we assume a trimer density of 10"^A"^. Moreover, since 
the effective intertrimer coupling sets the scale of the ordering 
temperature Tc, it is possible to reach values of Tc comparable 
or even higher than ambient temperature. 

Quantum magnets comprising weakly-coupled dimers be- 
came the focus of intense research during the last two 
decades | 37 , 38 1. We hope that the phase diagram presented in 
Fig. |2] will trigger a similar effort in the search for realizations 
of weakly coupled trimers. We note that typical perturbations 
of triangular molecules, such as deviations from equilateral 
shape or intratrimer Dzyaloshinskii-Moriya interactions, will 
not induce important changes in our phase diagram as long as 
they are small in comparison to f and f\ 
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No. DE-AC52-06NA25396 through the LDRD program. 



[1] C. Lacroix, R Mendels, and F. Mila, eds., Introduction to Frus- 
trated Magnetism: Materials, Experiments, Theory (Springer- 
Verlag, Berlin Heidelberg, 2011). 

[2] L. N. Bulaevskii, C. D. Batista, M. V. Mostovoy, and D. L 
Khomskii, Phys. Rev. B 78, 024402 (2008). 



4 



[3] K. L Kugel' and D. 1. Khomskii, Sov. Phys. Usp. 25, 231 
(1982). 

[4] G.-W. Chern and C. D. Batista, Phys. Rev. Lett. 107, 186403 
(2011). 

[5] H. U. Gudel, J. Chem. Phys. 82, 2510 (1985). 
[6] M. Honda, M. Morita, and M. Date, J. Phys. Soc. Jpn. 61, 3773 
(1992). 

[7] J. Padilla, D. Gatteschi, and P. Chaudhuri, Inorg. Chim. Acta 

260,217 (1997). 
[8] H. Sakurai, K. Yoshimura, K. Kosuge, F. Ishikawa, H. Mita- 

mura, and T. Goto, J. Phys. Soc. Jpn. 71, 664 (2002). 
[9] H. Lopez-Sandoval, R. Contreras, A. Escuer, R. Vicente, 

S. Bernes, H. Noth, G. JefFery Leigh, and N. Barba-Behrens, 

J. Chem. Soc, Dalton Trans, pp. 2648-2653 (2002). 
[10] B. Cage, F. Cotton, N. S. Dalai, E. A. Hillard, B. Rakvin, and 

C. M. Ramsey, C. R. Chim. 6, 39 (2003). 

[11] Y. Nakano, T. Yagyu, T. Hirayama, A. Ito, and K. Tanaka, Poly- 
hedron 24, 2141 (2005). 
[12] M. B. Stone, F. Fernandez-Alonso, D. T. Adroja, N. S. Dalai, 

D. Villagran, F. A. Cotton, and S. E. Nagler, Phys. Rev. B 75, 
214427 (2007). 

[13] K.-Y. Choi, N. S. Dalai, A. R Reyes, R L. Kuhns, Y. H. Mat- 
suda, H. Nojiri, S. S. Mai, and U. Kortz, Phys. Rev. B 77, 
024406 (2008). 

[14] J. Schnack, H. Nojiri, R Kogerler, G. J. T. Cooper, and 

L. Cronin, Phys. Rev. B 70, 174420 (2004). 
[15] N. B. Ivanov, J. Schnack, R. Schnalle, J. Richter, P. Kogerler, 

G. N. Newton, L. Cronin, Y. Oshima, and H. Nojiri, Phys. Rev. 

Lett. 105, 037206 (2010). 
[16] H. Manaka, T. Etoh, Y. Honda, N. Iwashita, K. Ogata, N. Ter- 

ada, T. Hisamatsu, M. Ito, Y. Narumi, A. Kondo, et al., J. Phys. 

Soc. Jpn. 80, 084714(2011). 
[17] M. Azuma, T. Odaka, M. Takano, D. A. Vander Griend, 

K. R. Poeppelmeier, Y. Narumi, K. Kindo, Y. Mizuno, and 

S. Maekawa, Phys. Rev. B 62, R3588 (2000). 
[18] Y. Qiu, C. Broholm, S. Ishiwata, M. Azuma, M. Takano, R. Be- 



wley, and W. J. L. Buyers, Phys. Rev. B 71, 214439 (2005). 
[19] H.-T. Wang, Phys. Rev. B 65, 024426 (2001). 
[20] S. Wessel and S. Haas, Phys. Rev. B 63, 140403 (2001). 
[21] S. Nishimoto and M. Arikawa, Phys. Rev. B 78, 054421 (2008). 
[22] K. Okunishi, S. Yoshikawa, T. Sakai, and S. Miyashita, 

Progress of Theoretical Physics Supplement 159, 297 (2005). 
[23] T. Sakai, M. Sato, K. Okunishi, Y. Otsuka, K. Okamoto, and 

C. Itoi, Phys. Rev. B 78, 184415 (2008). 
[24] M. Sato, Phys. Rev. B 75, 174407 (2007). 

[25] M. Sato and T. Sakai, Phys. Rev. B 75, 014411 (2007). 

[26] K. Okunishi, M. Sato, T. Sakai, K. Okamoto, and C. Itoi, 

arXiv: 1109.0063 (to be pubHshed in Phys. Rev. B). 
[27] M. Takahashi, P. Turek, Y Nakazawa, M. Tamura, K. Nozawa, 

D. Shiomi, M. Ishikawa, and M. Kinoshita, Phys. Rev. Lett. 67, 
746(1991). 

[28] Y. Nakazawa, M. Tamura, N. Shirakawa, D. Shiomi, M. Taka- 
hashi, M. Kinoshita, and M. Ishikawa, Phys. Rev. B 46, 8906 
(1992). 

[29] K. Shimizu, T. Gotohda, T. Matsushita, N. Wada, W. Fujita, 

K. Awaga, Y Saiga, and D. S. Hirashima, Phys. Rev. B 74, 

172413 (2006). 
[30] V. Subrahmanyam, Phys. Rev. B 52, 1133 (1995). 
[31] F. Mila, Phys. Rev. Lett. 81, 2356 (1998). 
[32] C. Raghu, I. Rudra, S. Ramasesha, and D. Sen, Phys. Rev. B 

62, 9484 (2000). 
[33] M. E. Zhitomirsky, Phys. Rev. B 71, 214413 (2005). 
[34] Here //sat is the field required to reach the 5^ = 1/2 state in 

every trimer. 
[35] Y. Kamiya and C. D. Batista, (unpubhshed). 
[36] J. H. R Colpa, Physica A 93, 327 (1978). 
[37] S. E. Sebastian, N. Harrison, C. D. Batista, L. Balicas, 

M. Jaime, P. A. Sharma, N. Kawashima, and I. R. Fisher, Nature 

(London) 441, 617 (2006). 
[38] T. Giamarchi, C. Ruegg, and O. Tchernyshyov, Nature Phys. 4, 

198 (2008). 



5 



